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A NUMERICAL METHOD FOR EVALUATING WAVE DRAG 
By Maurice S. Calm and Walter B. 01st ad 

SUMMARY 


A numerical method for evaluating the Von K&cm&n. wave-drag equation 
has heen developed and applied to the calculation of wave drag for sev- 
eral Bodies of revolution. Results indicated good agreement with the 
exact solution. Sxofficient accuracy of wave drag was obtained by using 
a simple numerical method to determine the second derivatives of the 
area distributions. 

It is concluded that the n-umerical method will yield results well 
within the accuracy of linearized theory. The method may be set up 
easily for a desk calculator or an electronic con^juter. 


INTRODUCTION 


Area-rule concepts (ref. l) have shown that the wave drag of a 
configuration is related to the wave drag of an equivalent body of 
revolution. As a result, much interest has been directed toward the 
evaluation of the wave drag of bodies of revolution. The most cranmon 
method of approach has been to eveiluate the Von K^rm&i wave-drag for- 
mula with a Fourier series analysis. This method is outlined in refer- 
ence 2. It woiild seem that a method of numerically evaluating the double 
integral in the Von KArmdn equation might, in some cases, be more useful 
to the engineer. Such a method is devised and presented in this report. 


SYMBOLS 


Cd 


wave-drag coefficient. 


Wave drag 
qS^ 


D wave drag 


i 


index of simimatlon in x 
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j index of svmmiatlon in x - | 

Lj defined by equation (k) 

I body length 

M free -stream Mach number 

n number of terms of sxjmmation 

q dynamic pressure 

r body radius 

E body maximum radius of configuration 1 

S body cross-sectional area 

Sf body frontal area of configuration 1 

X coordinate of longitudinal axis of body 

p = - 1 ’’ 

I auxiliary coordinate of longitudinal axis of body ^ 

Primes indicate derivatives with respect to the argument. 

ANALYSIS 

The Von Kdrmdn wave -drag equation for a body of revolution as given 
in reference 5 can be presented in the form 

1= - J S”(x)S"(|)log(x - g)d 5 dx (1) 

The integral in equation (l) may be considered as the volume between 
a surface determined by the function S”(x)S"(|)log(x - |) and the 
x,^ plane. The volume is bovinded laterally by the planes 6=0; x = I, 
and X = 6, as shown in the following sketch: 

u 
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Along any line x - | = Constant, the term log(x - |) Is a constant. 
Thus, if the integration proceeds first along this line, the term 
log(x - l) may be taken outside of the integral sign. The second 
integration is then performed with respect to (x - |) from 0 to 1. 

From these considerations, a ntaaerical solution to equation (l) can 
be somewhat simplified. The x, | plane can be divided into a number of 
finite diagonal strips of equal width, eind values of S"(x)S"(|) along 
the center of each strip can be con^juted. These values then are summed 
along the strips for which x - | = Constant and multiplied by the 
value of log(x - |) integrated across the strip. Using the integrated 
value of log(x - |) over the strip rather than the value of the loga- 
rithm itself' avoids the problon. of the singularity on the line x = | . 

It should be noted here that the summation for the line x = g is 
divided by 2 so that no areas outside of the limits of the integration 
are included. Finally, the products of the Btnnmation along each line 
for which x - | = Constant and the integrated value of log(x - g) 
are summed to obtain the solution. This integration is thus described 
by the following expression; 



n-1 n-1 


J=0 i=j 


s^s 


i-J 


( 2 ) 
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where Si' = S"r.i.fi + 

^ |n\ 2/ 

i^-¥ 


L4 = s 

J i 


for the exeitiples herein, and 


log(x - |)d(x - 1 ) 


- [(1 * * 5) - (■* - (■> - I)] 


> ( 5 ) 


Since S’(o) = S'(l) = 0, 
I r>x 


S"(x)S"(5)dg dx =|[S-(Z) - S'(0)] = 0 

Therefore the constant term ^log i in Lj can he eliminated, with 


the result that 


Lj = (j . i)los(j . I) - (j - |)log(j - i) 


(4) 


When J = 0, Lq = ^ log 


Equation (4) is Independent of both the number and the size of 
increments and can be used whenever S'(o) = S'(l) = 0. Values of the 
function Lj for j from 0 to 99 a^re presented in table I. 

When S’(l) is not equal to zero, the term ^log ^ must be 
retained in Lj, and additional terms must be used with Von Kdrmdn's 
equation ( see ref . 4) . These terms are 



log 


2 

Mil 


S'{1) 

jr 



S"(x)log(l 


- x)dx 


( 5 ) 


The integral in equation ( 5 ) can be evaluated with a single numerical 
summation by utilizing the information already obtained in the evalua- 
tion of equation (l) . 
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DISCUSSION 


In order to determine the accuracy of this numerical method, the 
wave drag of an analytical hody of revolution was computed hy this 

method, with the --and -i- - axes each divided into 40 equal Increments, 

II 

and by analytic integration of Von Karman’s equation. The shape of 
the analytical body (configuration l) was given by the following 
expression: 


^=4 


(!) - it 


In order to test the accuracy of the numerical procedure, exact values 
of the second derivative of the area distribution S" were used. The 
value of the wave -drag coefficient obtained by the numerical method was 

42.648^yj as compared with the exact value of 42.667^^^^, a difference 

of 0.045 percent. A layout of the calculations involved in the numer- 
ical method is presented in table II. 


The wave drag for configuration 1 was also determined numerically 
by using 100 increments. Again, exact values of S" were used. The 
value of the wave -drag coefficient obtained from these calculations was 


42.523(f)", 


a difference of 0.357 percent from the exact value. 


In practical applications of the numerical method, the exact values 
of the second derivative of the area distributions would not be avail- 
able. In fact, the exact area distribution is not generally known. 

Thus, evaluation of the second derivative by various numerical proce- 
diares may lead to considerable error. These errors, in turn, may have 
a large effect on the accuracy of the numerical method for evaluating 
wave drag. In order to determine this effect, three additional bodies 
of revolution were developed for which the exact values of S" were 
known. These bodies have large variations of curvature of their area 
distributions in order to provide a severe test. Configuration 2 was 


for 


0.4 ^ I < 


0.6 to 


obtained by adding 0.25 |l + cos 40jt^Y - 

the nondimens ional area distribution of configuration 1 (the parabolic 
body of revolution described previously) . Configuration 3 was obtained 
by subtracting this term from the nondimenslonal area distribution of 
configuration 1. Configuration 4 was obtained by adding 


0.25 


1 + cos - 0.4^ 


for 0.3 


^ T ^ 

L 


to the nondimenslonal 
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area distribution of configuration 1 . These area distributions (see 
fig. 1 ) were plotted to a scale conmensurate with the accuracy of the 
area distribution of a typical wind-tunnel model. Values of the second 
derivatives were then obtained by picking values of area from the curves 
and substituting them into the formula 



%-I - + Sn.1 

U/n)2 


where n = ICX). A comparison of these approximate values of S" with 
exact values is presented in figure 2 for configuration 5 - Values of 
wave-drag coefficient were computed for the four bodies by the numericEil 
method (for 100 increments) by using first, the exact values of S", 
and second, the approximate values. A ccraparison of the resTilts is 
shown in the following table: 



Cp 

for - 

Percent 

difference 

Configuration 

Exact 

S" 

Approxipiate 

S" 

1 

, /R \2 

42 . 525 (- 

\lj 

42.655(1)^ 

0.51 

2 

529.251^(1)^ 

525.815(1)^ 

1.04 

3 

505.255(1)^ 

296 . 81 t 8 (|)^ 

2.12 

k 

642.458(1)^ 

645. 914.9 

.55 


Despite the relatively large errors in some of the individual 
approximate values of S" (see fig. 2 ), the values of wave-drag coeffi- 
cient conputed from these values were in close agreement with those com- 
puted from the exact values of S". These results are not surprising 
when it is considered that the approximate values of S" form a bet of 
exact values of the second derivative of an area distribution which 
differs little from the original area distribution. The differences 
between the two area distributions will be of the same order of magni- 
tude as the accuracy to which the original area distribution is known. 
Obviously the wave drag determined by these two area distributions will 
be approximately the same. Any differences will be well within the 
aecinracy by which the theory can be expected to apply to a practical 
example . 
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The exaii5)les cited in the previous paragraphs indicate that 
100 increments are sufficient to yield good accuracy in the evaluation 
of wave drag. A larger number of increments might he used for an area 
distribution which has even more rapid changes in shape than those 
studied herein. However, it should be kept in mind that a body with 
such an area distribution will not permit linearized flow approximation 
and equation (l) should not be expected to yield good agreement with 
experiment. In fact, if the slope of the area distribution is discon- 
tinuous, Von Kfendti's equation Indicates an infinite value for the wave 
drag, which obviously disagrees with experimental evidence. 

It should be noted that the technique developed herein can be 
readily adapted to the evaluation of the wave drag of lifting conflgura 
tions (see ref. 5) said to vortex drag of a lifting surface in subsonic 
or supersonic flow (see ref. 14-). 


CONCLUDING REMAEKS 


A numerical method has been developed for evaluating Von K^rm^'s 
wave-drag equation. The method may be set up, easily for a desk calcu- 
lator or an electronic computer and will yield res\ilts well within the 
accuracy of linearized theory. A sinple numerical method was used for 
determining the second derivatives of nonanalytlc area distributions 
for foxar bodies of revolution. Resiilts of calculations made by using 
these approximate derivatives and by using exact derivatives yielded 
differences in wave drag on the order of 2 percent for a practical case 

The numerical method developed herein can be adapted to the evalua 
tion of the wave drag of lifting configurations and to the vortex drag 
of a lifting surface in subsonic or supersonic flow. 


Langley Aeronautical Laboratory, 

National Advisory Co 4 imlttee far Aeronautics, 
Langley Field, Va., February 28 , 1958 . 
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TABLE I.- VALUES OF THE FUNCTION Lj 
FOR i FRCM 0 TO 99 


J 


J 


j 


0 

-0.3^70 

34 

4.5265 

68 

5.2191 

1 

. 95^+7 

55 

4.5552 

69 

5.2346 

2 

1.6825 

56 

4.5855 

70 

5.2481 

5 

2.0939 

57 

4.6109 

71 

5.2631 

It 

2.5837 

58 

4.6377 

72 

5.2767 

5 

2.6078 

59 

4.6635 

73 

5.2898 

6 

2.7906 

40 

4.6888 

74 

5.3045 

7 

2.9^4-51 

4 l 

4.7155 

75 

5.3172 

8 

5.0788 

42 

4.7576 

76 

5.3309 

9 

3.1967 

43 

4.7613 

77 

5.3440 

10 

3.3022 

44 

4.7843 

78 

5.3567 

11 

5.3975 

45 

4.8065 

79 

5.5696 

12 

3.4846 

46 

4.8286 

80 

5.5820 

15 

3.5648 

47 

4.8505 

81 

5.5940 

14 

3.6387 

48 

4.8710 

82 

5.4071 

15 

3.7080 

49 

4.8919 

85 

5.4190 

16 

3.7724 

50 

4.9120 

84 

5.4304 

IT 

3.8331 

51 

4.9519 

85 

5.4431 

18 

3.8903 

52 

4.9512 

86 

5.4537 

19 

3.9443 

53 

4.9704 

87 

5.4664 

20 

5.9956 

54 

4.9^ 

88 

5.4770 

21 

4.0445 

55 

5.0072 

89 

5.4890 

22 

4.0909 

56 

5-0255 

90 

5.4997 

23 

4.1355 

57 

5.0428 

91 

5.5109 

24 

4.1779 

58 

5.0609 

92 

5.5218 

25 

4.2189 

59 

5.0776 

95 

5.5325 

26 

4.2579 

60 

5-0959 

94 

5.5455 

27 

4.2958 

61 

5.1113 

95 

5.5542 

28 

4.5322 

62 

5.1272 

96 

5.5637 

29 

4.3674 

63 

5.1429 

97 

5.5748 

30 

4.4010 

64 

5.1585 

98 

5.5855 

31 

4.4342 

65 

5.1747 

99 

5.5950 

32 

4.4656 

66 

5.1895 



33 

4.4963 

67 

5.2050 
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BLE n. - SAMPLE CALCULATIONS FOR PARABOLIC DODT OF REVOLUTION 

[la nuUdn; the calculations values of S" with seven decimal 
places were used, but for coavealeace of tabolatloa thoy have 
been rounded to four places] 





Bi 

B9 

E9i 

Oil 


Di 

IQ; 

^9i 

m 

m. 

ig 

^31 







1^ 


ji^ 



.2121 

•1793 

.1^ 

.OSOB 

.0661 

.0420 

.0197 

.0009 

.0196 

.0570 

.0525 

• .2Sp 

-.U506 

-.1139 

-.0634 

-.0330 

-.0246 

.0002 

.0230 

.0467 

.0662 

.0636 

-.3165 

-.2676 

-.2213 

-.3.779 

-.IW 

-.0963 

..062t 
-.0294 
.0014 
.0296 
• CW2 
.0782 
.0967 
.1167 

-.5579 

-.3006 

-.2507 

-.2014 

-.1530 

-.ms 

-.0709 

-.0332 

.0016 

.0335 

.0623 

.0663 

.mr 

.1320 

.1494 

-.3927 

-.3522 

-.2750 

-.2209 

-.1700 

-.1223 

-.0778 

-.0364 

,0017 

.0367 

.0663 

.09n 

.1446 

-.4204 

-.5337 

-.2944 

-.2565 

-.1621 

-.1310 

-.0655 

-.0390 

.0016 

.0393 

.0754 

.1040 

.1313 

.a»i 

.1733 

.19^ 

.2062 

■.3090 

-.2485 

-.1911 

-.1375 

-,067V 

-.0410 

.0019 

.0413 

,0770 

.1092 

.UT« 

.162s 

.1SW2 

.2021 

•21M 

.2271 

-.4552 

-.3^1 

-.3167 

-.2361 

-.1971 

-.1416 

-.0902 

-.(A22 

.0020 

.0426 

.<yf94 

.1126 

.1421 

.1679 

.1900 

.2063 

.2232 

.2343 

.2416 

-.4621 

-.3910 

-.5236 

-.2600 

-.2001 

-.1439 

.0020 

.0432 

.0606 

.1145 

.1445 

.1703 

.1909 

.2116 

.2266 

.2378 

.2455 

.2491 

i 

1 

-.4552 

-.3831 

-.3187 

-.2361 

-.1971 

-.1416 

.0020 

.0426 

.0794 

.1126 

.140 

.1679 

.1900 

.2003 

.2416 

.2453 

.2435 

.2416 

i 


I 
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11 


ft 



n-1 

b 

0-1 
h ^ 

l-J 


i-ev 

i-ej 


i-CT 

1.C8 

1-C9 

1-30 



1.^3 

1-04 


1-36 

1-3T 

1-30 


1 

l-J 


-.H2VL 

-.3866 

-.5V16 

-.2591 

-.2291 

-.1618 

-.0666 

-jxki 

.0659 

U834 

.2884 

.k 009 

.5339 

.6484 

.T83V 

-9E59 

::§S 

-.2565 

>.l£21 

-.3J10 

-^33 

-.0550 

.0018 

.0«3 

.073^ 

.loto 

.1513 

.i»i 

.1759 

.1968 

.2082 

.218V 

.2212 

.2266 

.2266 

.2232 

.2l6V 

.2062 

-.5927 

-.5522 

-.2750 

-.2209 

-.1700 

-.1223 

.0017 

.0367 

.0695 

.0971 

.1226 

.1^ 

.1639 

.179S 

.19^ 

.2Q£L 

.2005 

-gvtfi 

.2116 

.2065 

.2021 

.1926 

.1796 

-.5579 

-.30^ 

-.2507 

-.201V 

-.1550 

-.1115 

-.0709 

-.0332 

.0016 

.0555 

.0625 

.0685 

.1117 

.1320 

.1V9V 

.1839 

.1755 

.iav2 

.1900 

.1929 

.1929 

.1900 

U8V2 

.11.9* 

-.3163 

-.2676 

—2215 

--im 

-Jio9 

-0905 

-.0627 

-.029V 

.OOlV 

.0296 

.0552 

.0762 

.0907 

.1167 

.1320 

M6 

.1551 

.1626 

.1679 

.1705 

.1735 

.1679 

.1^ 

.1551 

.1446 

.1320 

.1167 

-.2677 

-.2265 

—IS^ 

-.1906 

-.1159 

-.0854 

-.0530 

.0012 

.0250 

.0467 

.0662 

.0636 

.0967 

.UX7 

A513 

.1376 

.1421 

.1443 

.1443 

.1421 

.1378 

.1513 

.1226 

.1117 

.0967 

.0636 

-.2121 

-.1795 

-.1485 

-n» 

-.0918 

-.0661. 

-.0420 

-.0197 

.0009 

.0196 

.0370 

.0525 

.0662 

.0782 

.0^ 

.0971 

U040 

.1092 

.1126 

.1143 

,1143 

.1126 

.1092 

.1040 

.0971 

.0685 

.0782 

.0662 

.0525 

-.1496 

-.1266 

-.1048 

-.0642 

-.064& 

-.04^ 

-.0296 

-.0139 

.0007 

J3140 

.0261 

.0370 

.0^7 

.0552 

.0683 

.0734 

.0770 

.07^ 

.0606 

J3Bo6 

.CFT94 

.0773 

.0754 

.0^ 

.0625 

.0552 

.0570 

.0261 

-.0347 

-.oe50i 

-.0139 

-.007V 

.0004 

.0075 

.0140 

.0196 

.0250 

.0296 

.0535 

.0367 

.0393 

.CVI3 

.0V26 

.0432 

.0432 

.0426 

.0413 

.0393 

.0367 

.0335 

.0^ 

.0250 

.QIS6 

.Q14C 

.0075 

-.0038 

-.0032 

-.0026 

-.0021 

-.0016 

-.0CE12 

-.0007 

-.0003 

0 

.0004 

.0007 

.0009 

.0022 

.0014 

.0016 

.0017 

.Ode 

.Od9 

.0020 

.0020 

.0020 

.0019 

.0018 

.0017 

.0016 

.0014 

.0012 

.0009 

.0007 

.oocv 

0 

.0798 

.0673 

.0557 

.0446 

.€045 

.0246 

.0156 

,00^ 

-.0003 

-.0074 

-UJ159 

-.01^ 

-.0246 

-.0294 

-.0332 

-.0364 

-.0390 

-.0410 

-.0422 

-.0429 

-.0429 

-.0422 

-.0410 

-.0390 

-.0364 

-.0532 

-.0294 

-.0248 

-.0197 

-.0139 

-.0074 

-.0003 

.0074 

.1699 

.1437 

.UB9 

.0956 

.0735 

.0529 

.0336 

.0156 

-.0007 

-.0159 

-.02^ 

-.0420 

-.o») 

-.0^ 

-.0709 

-.0778 

-.0833 

-.0074 

-.0902 

-.0315 

-.0915 

-.0902 

-.0674 

-.0633 

-.0778 

-.0709 

-.0627 

—0590 

-.0420 

-.0296 

-.0159 

-.0007 

.0058 

.0336 

.2671 

.2260 

.1673 

.1136 

.0632 

.0^ 

.0248 

-.0012 

-.0250 

-.0466 

-.0661 

-.0634 

-.0965 

-.1115 

-.1223 

-.1310 

-.:W5 

-.UklS 

-.1459 

-.1439 

-.1418 

-J375 

-.1310 

-.1223 

-.1115 

-.0995 

-0894 

-.0661 

-.0466 

-.0250 

-.0012 

.0246 

.0529 

.0032 

.3712 

.3141 

.2600 

.2069 

.1606 

.1156 

.0735 

.<045 

-.0016 

-.0547 

-.0648 

-.0918 

-ai59 

-.1369 

-.1550 

-.1700 

-.1621 

-.1501 

—1971 

-.2001 

-.2X1 

-.1971 

-.I9U 

-U821 

-.1700 

-.1550 

-0369 

-.1159 

-j39!16 

-.0646 

-.0347 

-.003^ 

.0345 

.0733 

.1156 

u6o6 

.46^ 

.4061 

.3578 

.2714 

.2069 

.1503 

-.0021 

-.0451 

-.0642 

-.1193 

-.1506 

-1779 

-.2014 

-.2209 

-.2463 

-.2561 

-.2600 

-.2600 

-.2561 

-.A85 

-.2365 

-.2209 

..2ca4 

-J7T9 

—1306 

-uss 

—oMa 

-j»5i 

-.0021 

.0446 

.0956 

.1505 

.2069 

.2714 

.6004 

.5ceo 

.4205 

.5578 

.2600 

.isro 

.1169 

.0057 

-.0026 

-.0561 

-.1046 

-.1485 

-,l8T-» 

-.2215 

-fi507 

-.2750 

..2944 

-.3090 

—5187 

-.5236 

-.3236 

-.3107 

-.3090 

-.2944 

-.2750 
-.2507 
-.ai5 
— iflTV 
-.1485 
-.1048 
-.C56I 

-W0026 

.0557 

.1189 

.1070 

.2600 

.3578 

.4205 

-725V 

.6138 

.5060 

.4061 

.5141 

.2260 

.1437 

^75 

-.0032 

-.0676 

-.1266 

-.2676 

-.3026 

-.3322 

-.5557 

-5754 

-.3051 

-.5910 

-.3910 

-.3051 

-.3734 

-.3537 

-.3322 

-.3026 

-’.1795 

-.1266 

-.0678 

-.0032 

.0673 

.1437 

.ss6o 

.5141 

.4061 

.5060 

.6138 

.0574 

.7254 

.6004 

.3712 

.2671 

.1699 

.0796 

-.0038 

-.0002 

-.1496 

-.2121 

-.3579 

-.3927 

-.4204 

-.4413 

-.4332 

-.4621 

-.4621 

-.4552 

-.4413 

-.4204 

—3927 

-.3579 

-.5163 

-.2677 

-.2131 

-.1496 

-.C602 

-.0038 

.0^ 

.16» 

.2671 

.3712 

.4624 
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Figure 2.- Comparison of exact and. approximate values of the second derivative of the area dis- 
tribution of configuration J. The solid line indicates the exact values of S". The symbols 
indicate approximate values of S" . 



